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Collapse models postulate the existence of intrinsic noise which modifies quantum mechanics and is re¬ 
sponsible for the emergence of macroscopic classicality. Assessing the validity of these models is extremely 
challenging because it is non-trivial to discriminate unambiguously their presence in experiments where other 
hardly controllable sources of noise compete to the overall decoherence. Here we provide a simple procedure 
able to probe the hypothetical presence of the collapse noise with a levitated nanosphere in a Fabry-Perot cavity. 
We show that the stationary state of the system is particularly sensitive, under specific experimental conditions, 
to the interplay between the trapping frequency, the cavity size, and the momentum diffusion induced by the 
collapse models, allowing to detect them even in the presence of standard environmental noises. 
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I. INTRODUCTION 


Quantum mechanics (QM) has been proven as an extremely 
accurate theory for describing objects at microscopic scales. 
It intrinsically shows no limit to describe large massive sys¬ 
tems GM- However, the lack of observations of macroscop- 
ically distinguishable superposition states of macroscopic ob¬ 
jects allows to conjecture that QM should be modified at large 
scales. The modifications are aimed to explain the collapse of 
the wave function at the macroscopic level and also to solve 
the quantum measurement problem Q. Various attempts 
have been made in this direction, including the Ghirardi- 
Rimini-Weber approach mi , continuous spontaneous localiza¬ 
tion (CSL) 121, and gravitationally-induced collapse models 
(CMs) ISl. These models modify the Schrodinger equation by 
introducing appropriate stochastic non-linear terms, which re¬ 
solve the problems at macroscopic scales, while reproducing 
the standard results at microscopic scales E). 

Whether the proposed CMs are or are not exact should be 
examined by experiments. Matter-wave interferometry oni, 
where large massive molecules or clusters IfTTI are sent 
through interference gratings, is considered as one of the ideal 
arenas to test CMs. However, the mass range has, so far, yet to 
be reached to effectively test CMs. An alternative approach is 
based on cavity optomechanics 112, where one could prepare 
massive mechanical resonators (MRs) in linear superposition 
states and monitor their decoherence lEHBl- More recent 
proposals lfT6H20ll suggested to test CMs in optomechanical 
systems in a non-interferometric way, so that the preparation 
of large spatial superposition states is not required. In fact, 
the spontaneous collapse mechanism leads to spatial deco¬ 
herence, i.e., momentum diffusion, of the MR, which results 
in additional phase noise of the light leaking out of the cav¬ 
ity OMISl. 

In general, the identification of systems and regimes in 
which collapse-induced diffusion is theoretically dominant 
over the environmental noises is not sufficient for the design 


of experiments able to univocally decide whether a given ob¬ 
servation is actually the result of collapse theories, or of other 
uncontrolled sources of environmental decoherence. Being 
able to differentiate unambiguously their effect is a non-trivial 
task which deserves further studies and the identification of 
dedicated experimental procedures. The fundamental idea at 
the base of this work, is the observation that different sources 
of noise exhibit different scalings with the system parameters, 
and hence distinguishable scalings of measurable quantities 
are expected to be observed if determined sources of noise are 
or are not actually present or effective. Specifically, in this ar¬ 
ticle, we discuss a novel and efficient test of the CSL model 
in optomechanical systems in the regime of high mechanical 
quality factors and cryogenic temperatures, that is realizable 
with trapped levitated nanospheres in Fabry-Perot optical cav¬ 
ities maiiiMi. We demonstrate that, in these systems, the 
different noise sources are particularly sensitive to the trap¬ 
ping frequency co and to the length of the cavity L, so that the 
validity of the CSL model can be actually probed by the study 
of the nanoparticle dynamics as a function of co and L. This 
observation is general and can be applied to any optomechan¬ 
ical scheme. In the following, we explore its effectiveness in 
the analysis of the system steady state, by analyzing the phase 
noise of the light field leaking out of the cavity as a function 
of o) and L. We show that, in experimentally achievable pa¬ 
rameter regimes, a nanosphere can be prepared in a stationary 
state particularly sensitive to the mechanical momentum dif¬ 
fusion induced by the CMs, and, most importantly, whose sta¬ 
tistical properties scale differently depending on whether CMs 
are true or not, hence allowing for an efficient test of CMs. 


The remainder of this paper is organized as follows: In 
Sec.j^we describe in detail our model, provide the Langevin 
equations governing the dynamics, and analyze various con¬ 
ventional diffusion rates of the system as well as the nontrivial 
diffusion rate induced by the collapse noise postulated in the 
CSL model. In Sec. IIA and |IIB| we provide two alternative 
proposals to effectively test the CSL model, i.e., by varying 
the trapping frequency and the cavity length. We show de- 
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tails on the pressure and temperature, as well as the degree 
of precision for the measurement, required for testing the cor¬ 
responding value of the collapse rate. Finally, in Sec. Ill we 


draw our conclusions and in the Appendix we discuss the ef¬ 
fects of blackbody radiation in our system. 


II. THE MODEL 


We consider a single nanosphere of radius R trapped by a 
harmonic dipole trap, at frequency to, within a Fabry-Perot op¬ 
tical cavity with length L, finesse 'F, and with mirror radius of 
curvature R^. A single cavity mode, with resonance frequency 
ojc, is driven by an external field at power P and frequency 
a>L, detuned by A = cuc - cul, and is coupled to the center of 
mass of the nanosphere uniEBEa. The relevant degrees of 
freedom for the linearized system dynamics are the fluctua¬ 
tions of the cavity field and of the mechanical center of mass 
variables about their respective average values, described by 
the bosonic operators a and (with [a, a^] = 1) for the cavity 
field, and by the dimensionless position and momentum x and 
p (with {x,p} - i) for the nanoparticle. The corresponding 
quantum Langevin equations (QLEs) read 

a - - + k) a - \GX + , (1) 

X - cop, 

p — —ux-yp — G{a + a'^ + Fau + FD, 

where the linearized coupling strength G = g a is proportional 
to the average cavity field a - with 

g the bare optomechanical coupling, which can be expressed 

im, with Ac the cavity wavelength, 
e the electric permittivity of the nanosphere. Vs its volume, 
and Vc-ttLWqIA the cavity mode volume with mode waist 
VTo = [AcL(2Rc/. K-ncKl^FL) is the cavity 
linewidth with c the speed of light, and y is the damping 
rate of the mechanical motion. For levitated nanospheres y 
can be extremely small, resulting in very high quality factors, 
> 10^° ll26l . with the dominant contribution due to friction 
from residual air molecules, for which y - with Pa 

the gas pressure, v - ^IJkdTjnTa the mean speed of the air 
molecules, lUa their mass (which we take nia - 28.97 amu), 
and T the air temperature 11211 . In Eq. Q we have included 
the relevant sources of noise affecting the dynamics of the sys¬ 
tem, and leading to mechanical Brownian motion, in terms of 
the 5-correlated stochastic forces a‘"(f), Fair(f) and Foit) (see 
the Appendix for a comment on the effects of blackbody radi¬ 
ation). Eirstly, «”(?) is the input noise operator for the cavity 
field due to the fluctuations of the external electromagnetic 
environment. Its only non-zero correlation function is 

{d'^{t)d^\t'))^6{t-t'). ( 2 ) 

The term F^idt), instead, accounts for the mechanical noise 
due to the scattering of background air molecules which is 
related to the dissipation rate y by the fluctuation-dissipation 


theorem, and its autocorrelation function, in the relevant high 
temperature regime, is given by 

{F,iAt)F,isit'))^Da6{t-t') (3) 

with the diffusion rate Da = Fyk^Tl(Hto). Einally, Fo(f) ac¬ 
counts for pure diffusion of the particle motion with autocor¬ 
relation function 


(Foit) Foil')} = (£» + dsph)5(f - 0, (4) 


responsible for dephasing and decoherence on the nanosphere. 
Here we have separated the contributions D and Tsph which 
describe, respectively, the effects of light scattering and of 
collapse-induced diffusion. 

In the framework of the CSL model, the collapse-induced 
diffusion rate for a spherical nanoparticle, with constant mass 
density po and harmonically trapped at frequency co, is given 

byUni 


A pii 

'^sph “ T 

CO 


m: 


0 




2r? 


+ 1 ) 


(5) 

R^ 


with mo the atomic mass unit. The actual strength of collapse 
noise is determined by two phenomenological parameters, the 
characteristic length and the collapse rate A. While there is 
significant agreement on the estimated value for the charac¬ 
teristic length Kc - 100 nm, the expected value of A is more 
controversial. The initial estimate of A is 10“'® s ‘ EH, 
however, larger values have been proposed by other authors 
(for example 10“**^ s“* in Ref Il27ll '). In any case, differ¬ 
ent experimental results provide indications that A should be 
lower than 10“*' s“' ES), 10“® s“' 1291 and 10“" s“' l30l . 
for rc - 100 nm. Relevant new proposals should be able to 
confirm or improve such results by either lowering the upper 
bound of A or by detecting noise effects which cannot be ex¬ 
plained by standard decoherence. 

The diffusion rate via photon scattering, instead, can be ex¬ 
pressed as the sum of the contributions due to the scattering 
of trapping and cavity light, i.e., D-D,+Dc, which are given, 
respectively, by EH 

' 9poco oju ’ ' 9pow Vc ' ^ ’ 

where ec - 3^, kc - 2nlAc, cou is the frequency of the 
trapping laser, and the trapping frequency is determined by 
( 1 ) - {AecII(pQcWf)YI^, with VT, the waist of trapping light 
which can be approximated by Wt ^ AdiTcN) with N the 
numerical aperture, and I the intensity of the trapping held, 
which is related to the power by J = PtHjcWj). Einally Uph = 
\a\^ is the mean cavity photon number. 

We observe that D, and Dc increase while Da decreases 
with the size of the nanoparticle. We have checked that the op¬ 
timal size of the sphere for testing the CSL theory is roughly 
around HU. We also note that large held powers imply 
large diffusion rates by light scattering and as a result low 
powers are in general required to reduce the photon scatter¬ 
ing induced diffusion and to make it comparable to the col¬ 
lapse noise, entering therefore a regime in which CSL could 










3 


be detected. These observations are however not sufficient to 
discriminate the effect of collapse induced diffusion. In order 
to achieve this, we have to identify strategies which exhibit a 
qualitatively different response when CSL is present and when 
it is not. In this respect, the central observation of this work is 
that the diffusion rates exhibit very peculiar scalings as a func¬ 
tion of a> and L. As we will discuss below different scalings 
of the diffusion rates imply distinguishable steady state behav¬ 
ior of the nanosphere which may be exploited to distinguish 
the effects of CMs. Specifically, in our setup, where we use a 
weak driving field at a frequency smaller than the cavity fre¬ 
quency which yields a weak cooling force on the mechanical 
motion and stabilizes it, the position fluctuation gets encoded 
in the phase quadrature of the cavity field Y - -i(a - a^')/ V2, 
such that (jc^) oc (T^) -h shot noise. Consequently direct in¬ 
formation on the nanosphere diffusion can be extracted by the 
measurement of the optical phase at the cavity output. Here¬ 
after we analyze the steady state behavior of (Y^) in various 
experimental conditions, versus either w or L, which may be 
used to probe the strength of A. 


A. Test of the CSL model by varying the trapping frequency 


We first note that A^j,h,Da oc l/a>, D, oc cu, and, when G is 
fixed, Dc is independent from w, as clearly shown in Fig. [T] 
(a). Specifically, Dc can be made negligible for sufficiently 
small G, i.e., for sufficiently small driving power. Similarly, 
under the conditions of low pressure and temperature. Da is 
much smaller than Tsph and D,. This is the situation achieved 
in Fig. (b), where we consider a larger size of the sphere 
R - Tc, in order to further reduce Da and increase Tsph. Thus, 
if in addition the power of the trapping light is small enough 
(corresponding to relatively small w), so that D, and dsph are 
of comparable strength, then the presence of spontaneous col¬ 
lapse mechanism can be demonstrated by detecting the out¬ 
put light as a function of the mechanical trapping frequency, 
which could be spanned by simply adjusting the intensity of 
the trapping light. In detail, the steady state variance of the 
optical phase quadrature displays distinguishably different be¬ 
havior depending upon the presence or absence of collapses. 
This is shown in Fig. [^(d) and Fig. where we analyze the 
results for a diamond nanosphere with radius = 100 nm and 
a cavity with finesse T - 10^, and we compare the results 
with (blue lines) and without (red lines) the effect of CSL. In 
particular, the stationary variance (F^) in the presence of the 
CSL effect increases rapidly as (jj is gradually reduced, while 
it is practically independent upon w without the CSL effect. 
Therefore, by repeating the experiment at different trapping 
light intensity one could verify this different behavior and de¬ 
termine the possible presence of spontaneous collapses. When 
the value of A is reduced as in Fig. the relative difference 
between the two curves (with and without the CSL) reduces. 
At /I ~ 10“*^ s“', the two curves become hardly distinguish¬ 
able, implying that the experimental realization of our proto¬ 
col would allow, if no CSL effects are detected, to lower the 
upper bound of A to 10 s '. 






a> (units of k) 


to (units of k) 


FIG. 1: (a), (b): Diffusion rates for the scattering of trapping light 
D, (blue), cavity light (red) and air molecules Da (gray), and the 
collapse rate dsph (green) versus the mechanical frequency tu with (a) 
R = Q.Svc and G = 0.05*:; (b) R = Vc and G - O.Ol*. Note that 
the curves of Dc and Da in (b) are very close to the tu-axis and no 
longer visible, (c), (d): Steady state variance of the optical phase 
quadrature (F^) versus the trapping frequency tu. Blue (red) lines 
refer to the case with (without) the CSL effect. The parameters for 
(c)/(d) correspond to those for (a)/(h). The other parameters are L = 
Rc - 1 cm, T = 10^ (corresponding to *• = 0.47 Mffz), A = 0.01*:, 
Ac = 1064 nm, N = 0.6, T = 1 K, = 10-‘° Torr, A = lO^* s^', 
Vc = 100 nm, and we consider a diamond nanosphere with po = 3.5 
g/cm^ and e - 5.76. 


B. Test of the CSL model by changing the cavity length 

An alternative approach provides a more evident effect of 
collapse noise at the expense of a slightly more involved ex¬ 
perimental protocol. In particular, we find that the diffusion 
rates show peculiar scalings with the size of the cavity L when 
the ratios w/*:. A/*: and G/*: are kept fixed. Specifically, since 
*■ cx 1/L then Tsph oc L and D, oc 1 /L, while Dc oc i2RclL- 1 
decreases with L due to the particular scaling of the geometry 





FIG. 2: Steady state variance of the optical phase quadrature (F^) 
versus the trapping frequency cj for (a) A = 10“^ T = 200 mK, 
Pa = lO-'o Torr; (b) A = IQ-'® s-‘, T = 100 mK, Pa ^ 3 x lO^" 
Torr; (c) A = 10^" s^', T = 60 mK, Pa = 10^“ Torr; (d) A = lO^'^ 
s“', r = 10 mK, Pa = 10“'^ Torr. Blue (red) lines represent the case 
with (without) the CSL effect. In all plots we take G = O.OOl/c. The 
other parameters are as in Fig. □(b). 
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of the cavity mode. In this case the experimental procedures 
should run as follows. For each value of the cavity length 
(implying different values of the cavity linewidth) the values 
of u>, G and A should be carefully tuned and monitored in 
order to achieve determined fixed values relative to k, so that 
the system remains under the same optomechanical condition. 
This can be achieved by varying the intensity and frequency of 
the stabilized driving and trapping laser, and simultaneously 
monitoring the G and a> values. The resonance frequency a> 
is easily measured from the position of the resonance peak in 
the cavity output phase, while G can be extracted from the 
nanosphere cooling rate, which is given by G^/a: in the bad 
cavity regime considered here. The steady state observables, 
detected for each set of parameters, will thus depend only on 
the diffusion rates. Thereby, if the powers of both trapping and 
driving light are sufficiently small (corresponding to relatively 
small to and G), so that the D,, Dc and Tsph are of comparable 
strength, then the behavior of the system steady state versus L 
may distinguish the action of CMs. This is shown in Fig. 
where the steady state value of (T^) is reported versus the mir¬ 
ror distance, with (blue lines) and without (red lines) the effect 
of CSL. In the plots, the cavity length L - corresponds to a 
confocal cavity, whereas the largest value of L approaches the 
limit of a concentric cavity L - 2Rc. We remark that the cav¬ 
ity is unstable for larger L, so that the presented results cover 
all the possible geometric configurations of a stable symmet¬ 
ric Fabry-Perot resonator. We observe, in Figs. (a) and (b) 
that as the trapping frequency to and the linearized coupling G 
are gradually reduced, the effect of the collapse-induced dif¬ 
fusion becomes more and more distinguishable. Specifically, 
the presence of spontaneous collapses is signaled by a change 
in the slope of these curves for L > Re- When the value of A 
is decreased, as in Figs.[^(c) and (d), the effect of CSL be¬ 
comes less and less visible, although here results similar to 
those reported in Fig. can be achieved at a slightly larger 
temperature. In particular, at d = 10“'^ s“', 10“^' s“' and 
10“*° s“* the relative difference between the two curves, when 
L is large, is of ~1.5%, ~12% and ~30%, respectively. This 
implies that both the precision of measurement of the optical 
phase and the precision with which the experimental parame¬ 
ters are calibrated and kept fixed in the repeated experimental 
runs must be smaller than the above relative difference. That 
is, with a parameter calibration and measurement precision of 
no more than ~1.5%, one could discriminate the CSL down to 
A = 10“*2 S-*. 


III. CONCLUSIONS 

The presented results suggest a general strategy to probe the 
possible effect of the CSL model and to discriminate it from 
the effect of other sources of decoherence. It is designed to 
work with levitated nanospheres in optical Fabry-Perot cavi¬ 
ties. The very high mechanical quality factor of these systems 
makes them the ideal platforms to test the CMs. Moreover, 
our approach is effective when the pressure and temperature 
of background gas are sufficiently small in order to make the 
corresponding noise negligible. In particular, we have dis- 
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L (units of Rc) 


FIG. 3: Steady state variance of the optical phase quadrature (V^) 
versus the cavity length L with (a) w = A = 0.02a', G = 0.15a, = 

10-*“ Torr; (b) to = A = 0.008a, G = 0.08a, = 10-*“ Torr; (c) 

w = A = 0.008a, G = 0.025a, = lO^’* Torr; (d) to = A = 0.008a, 

G = 0.009a, Pa = 10“'^ Torr. The values of A are reported in each 
plot. All the curves are evaluated for T = 100 mK and = 2 cm, 
and the other parameters are as in Fig. Blue (red) lines represent 
the case with (without) the CSL effect. 


cussed results for the optical phase quadrature Y and we have 
demonstrated that the presence of the CSL effect can be tested 
by investigating the stationary behavior of (T^) as a function 
of the trapping frequency to and of the cavity size L. Firstly, 
we have shown that when G is sufficiently small, the results 
without CSL are basically independent from u, while they in¬ 
crease rapidly as to decreases when the collapse noise is taken 
into account, hence providing a distinct signature of the CSL 
effect. Then we have also shown that by tuning the cavity 
length and proportionally also the field powers and frequency, 
the ratios between the optomechanical parameters (frequen¬ 
cies, coupling and photon loss rate) can be kept fixed, hence 
keeping the system in the same optomechanical regime, while 
the diffusion rates by light scattering and by CSL exhibit di¬ 
verging behavior. As a consequence the results for (T^) reflect 
a similar behavior, and thus the effect of the CM can be clearly 
discriminated. 

In our analysis we have considered an optomechanical sys¬ 
tem comprising a diamond nanosphere of radius 100 nm 
trapped, by an optical dipole trap of a few kilohertz, inside a 
Fabry-Perot cavity with finesse of 10^ and length of a few cen¬ 
timeters. Similar systems have been described in Refs. Il22l - 
l25l . We have thereby demonstrated that these protocols can 
be employed to test the strength of A to values as low as 10“*^ 
s“* with realistic parameters. This value is essentially limited 
by the temperature and pressure that can be achieved in exper¬ 
iments (we have considered values as low as - 10“*^ Torr 
and T = 10 mK, in the results versus to, and T = 100 mK, in 
the results versus L) and lower values of T and Pa would allow 
to test even lower values of A. These values of temperature and 
pressure, although challenging, have been already discussed 
in various experiments with cold atoms ll^ and are expected 
to be achievable in the near future also in experiments involv¬ 
ing nanospheres. We have proved that these results can be 
obtained by the measurement of the phase quadrature of the 
cavity held, which can be accessed with standard optical tech- 
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niques. On the other hand, we remark that similar consid¬ 
erations and results are in principle valid for any observable 
of the optomechanical system. This unique versatility makes 
the presented proposal very promising for an actual test of the 
CSL theory. 
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APPENDIX 

In this Appendix we briefly discuss the effects of black- 
body radiation in our model. In our treatment we have ne¬ 


glected the effect of blackbody radiation which is very small 
in the parameter regime addressed in this paper. In general 
the trapping and cavity light can heat up the particle and the 
corresponding emitted black-body photons can act as a noise 
source for the particle motion (absorption of black-body pho¬ 
tons is less relevant) dEI]. This effect can be estimated 
using the approach described in the Supporting Information 
of Ref. im. Here we consider diamond which is transparent 
from ultra-violet to infrared wavelengths, with small absorp¬ 
tion due to the quantity and the quality of possible impurities. 
Using the theory of Ref. ||2TI| with a relatively large absorption 
coefficient of 1 m * at the cavity field wavelength and a com¬ 
plex relative permittivity with an imaginary part of 10“^ at the 
black-body wavelength, we find an internal temperature of the 
particle of a few hundred degrees Kelvin, and a corresponding 
diffusion rate that in the worst case is three order of magnitude 
smaller than the smallest diffusion rate that we have included 
in our description. 
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